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Abstract 

We compare the light hadron spectrum and decay constants for quenched QCD 
at (5 = 6.2 using an 0(a)-improved nearest-neighbour Wilson fermion action with 
those obtained using the standard Wilson fermion action on the same set of 18 gauge 
configurations. For pseudoscalar meson masses in the range 330-800 MeV, we find no 
significant difference between the results for the two actions. The scales obtained from 
the string tension and mesonic sector are consistent, but differ from that derived from 
baryon masses. The ratio of the pseudoscalar decay constant to the vector meson mass 
increases slowly with quark mass as observed experimentally. 
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1 Introduction 



The numerical simulation of QCD involves systematic errors arising from a variety of sources: 
non-zero lattice spacing, finite volume, quenching and/or unphysical quark masses. Here we 
address the first of these by working in the quenched approximation at a fixed volume, and 
comparing the results obtained on the same set of gauge configurations using two formula- 
tions of the fermion action with different discretisation errors. 

Classically, the standard Wilson pure gauge action differs from the continuum Yang-Mills 
action by terms of 0(a 2 ), where a is the lattice spacing [|lj, whereas the Wilson formulation 
of lattice fermions introduces an 0(a) discretisation error in order to avoid the doubling 
problem. At the quantum level, matrix elements computed with the Wilson action have 
errors of 0(a). Monte-Carlo determinations of the renormalisation constants of the vector 
current have shown that these discretisation errors may be as large as 30% 0. 

In the spirit of the Symanzik improvement programme ||, which sought to eliminate in a 
systematic way discretisation errors at the quantum level, Wetzel [Q proposed a two-link 
fermion action which cancels the 0(a) term in the Wilson action at tree level. In a study of 
on-shell improvement, Sheikholeslami and Wohlert || introduced a nearest-neighbour 0(a)- 
improved fermionic action, which is more convenient for large scale simulations and which 
we refer to as the 'SW action. This may be obtained from the two-link action by a rotation 
of fermion fields in the functional integral. Heatlie et al. |2|] demonstrated the absence of the 
leading logarithmic terms at n th order in perturbation theory, in matrix elements of rotated 
operators. These terms are of the form g^aln 71 a, and hence effectively 0(a) in the weak 
coupling regime, because g$ ~ 1/ In a. 

In this paper we present further details of a study in quenched QCD on a 24 3 x 48 lattice at 
(3 = 6.2, using both the standard r = 1 Wilson fermion action: 



X )Q ~ l^) U ^(x)q(x + p) + q(x + //)(! + ^UUx)q{x) 



x { n L 

and the SW fermion action ||: 

S S F W = - z~ £ q(x)F, u (x)a, u q(x). (2) 

* X,fJ,,U 

F^ v is a lattice definition of the field strength tensor, which we take to be the sum of the 
four untraced plaquettes in the \w plane open at the point x ||, as indicated in Figure [|: 



F^(x) = [U n , u (x) - U^ lu (x) 

For this reason, the SW action is sometimes called the 'clover' action 



(3) 



Our objective is to look for evidence of improvement in masses and decay constants of the 
light hadrons. A summary of our results has appeared in 0, where we concluded that the 




Figure 1: Lattice definition of the field strength tensor, F flll (x) 



SW action for medium to light quark masses gives results in agreement with the Wilson 
action, although typically with a noisier signal. In this paper, we report our attempts to 
enhance the signal when using the SW action by smearing the same set of propagators at 
the sink. 

The plan of the paper is as follows. We begin by describing the generation of the gauge 
fields, and the measurement of the static quark- ant iquark potential used to extract the 
lattice spacing in physical units from the string tension. We describe the algorithms used 
to compute the quark propagators, and the bootstrap method used to analyse the hadron 
correlators. We compare the meson and baryon spectra for the two actions, and extract 
values for the lattice spacing by comparing the masses to their physical values. We investigate 
the hyperfine mass splittings as quantities that might be sensitive to the difference between 
the actions. Finally, we compute the pseudoscalar and vector meson decay constants. 



2 Computational Procedure 

The gauge field configurations and quark propagators were obtained using the Meiko i860 
Computing Surface at Edinburgh. This is a MIMD system of 64 nodes, each of which consists 
of one Intel i860 processor on which the application runs, two Inmos T800's for inter-node 
communications and 16 MBytes of memory. The peak performance is 5 Gflops in 32-bit 
arithmetic. Our codes are written in ANSI C with key numerically-intensive routines in 
i860 assembler, and communications using Meiko CS-Tools J7L . Apart from certain global 
operations which use 64 bits, the calculations are performed in 32-bit arithmetic. We achieve 



a maximum performance of 3 Gflops for some routines, while overall our programs sustain 
around 1.5 Gflops ||. 



2.1 Gauge Configurations 
2.1.1 Update algorithm 

There are computational advantages in using an over-relaxation (OR) algorithm for updat- 
ing, both from the speed of the code itself and from the improved transport through the 
space of equilibrium configurations. For SU (2) pure gauge theory, such an algorithm is very 
simple and hence computationally fast. For SU(3) pure gauge theory, there is no such sim- 
ple implementation known. The proposal of Creutz |10[ uses an approximate over-relaxation 
which has to be corrected by an accept/reject calculation. Furthermore, a projection of a 
3x3 complex matrix to SU(3) is needed, which is computationally slow. 

As an alternative, we explored the possibility of using OR steps in a sequence of SU(2) 
subgroups of SU(3). This is in the spirit of the Cabibbo-Marinari |Tl| scheme of updating 
SU(3) via SU(2) subgroups. Our SU(2) subgroup OR scheme is exactly micro-canonical, 
unlike the Creutz scheme. In order to optimise this SU(2) subgroup OR approach, we mea- 
sured the average distance moved by a link in the SU(3) group manifold (along a geodesic) 
for different subgroup updating schemes. Using 3 subgroups (12, 23 then 31) was optimal 
and gave a distance moved slightly less than the Creutz scheme, but significantly more than 
the Cabibbo-Marinari heat-bath update. A check of auto-correlation between configurations 
was made and confirmed that this SU{2) OR scheme is competitive. Because of its compu- 
tational speed, we chose to use this 3 subgroup SU (2) OR algorithm, complemented by one 
heat-bath update every 5 OR updates to preserve ergodicity. We call this type of algorithm 
'Hybrid Over- Relaxed'. At the end of each 5 + 1 sweeps, we re-unitarise each element of the 
gauge configuration by normalising the first row of the SU(3) matrix, making the second 
row orthonormal to the first and then reconstructing the third row from the appropriate 
vector product. We have checked that any systematic effect arising from this procedure is 
negligible. 

Whenever the gauge configuration is written to disk, the last operation performed on it is 
a re-unitarisation. This allows us to write only the first two rows of each matrix and yet 
re-start from the identical configuration by reconstructing the final row. 



Random numbers were generated using the uni random number generator ||12|| . Each pro- 
cessor used a separate copy of the generator, initialised with different starting states. 

The results presented here are based on an analysis for each action of the same set of 18 
configurations, starting at configuration 16800 and separated by 2400 sweeps. 



2.1.2 Static quark— antiquark potential 



The potential between static colour sources at separation R can be obtained by measuring 
Wilson loops of size R x T. However, it is useful to consider generalised correlations of size 
RxT where the spatial paths are any gauge invariant path with cylindrical symmetry (Ai g ), 
not just the straight line used in a rectangular Wilson loop. Moreover, the two spatial paths, 
at times T apart, may be different. So the measurable is Cij(R,T) where i, j label the path 
type. Then the largest eigenvalue, X(R), of the transfer matrix in the presence of the static 
sources is related to the potential by 

V(R) = -log(X(R)) (4) 

and X(R) can be obtained from the limit as T — > oo of the largest eigenvalue, X(R,T), of 
the equation 

Cij(R, T)uj = X{R, T)Cij(R, T - l) Uj . (5) 

This is a variational method where an optimal combination of paths i (with i — 1, ..N) is 
chosen. The best choice of paths would allow X(R,T) to be close to its asymptotic value for 
small T. In this way the statistical errors are minimised. Thus it is convenient to introduce 
an 'overlap' defined as 

OW = H{\{R,T)/\{R)}. (6) 

T 

For R = 12 straight paths, we find O < 0.006. Thus rectangular Wilson loops are a very 
inaccurate way of extracting the ground state potential. 

As has been known for some time, an efficient improvement comes from fuzzing or blocking 
these paths [jOH , |T4j 1~5| . Using a purely spatial blocking allows the transfer matrix interpre- 
tation to be retained. Because of its flexibility, we used the recursive blocking scheme [[Uj: 



[/(new) = Vsu(3) cCZ(straight) + £/(u-bends) . (7) 

i 

After exploratory studies to optimise the overlap O, we chose c = 4 and a maximum of 40 
iterations of this recursive blocking. As a variational basis with N = 2, we used 40 and 28 
iterations for the paths. This basis for R = 12 yields an overlap O = 0.95 which is a huge 
improvement over the purely unblocked case. Because this value of the overlap is close to 
1.0, one can get reliable estimates of the ground state potential from quite small T-values 
such as T = 3. A powerful cross check of this extraction of the ground state energy at rather 
small T-values comes from determining the first excited state also in the variational method. 



From previous work, this first excited A\ g potential should lie about 2n/R higher [15] and 
this is completely consistent with our results. 

To estimate the potential V(R), we need to extrapolate in T. For 1 < R < 4, we find that 
the 3 : 2 and 4 : 3 T-ratio results for X(R,T) are consistent within errors. As a conservative 



R 


V(R) 


V(R+1)-V(R) 


1 


0.3775(2) 


0.1564(3) 
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0.5339(3) 


0.0787(5) 


3 


0.6126(7) 


0.0509(5) 


4 


0.6635(8) 


0.0416(12) 


5 


0.7043(17) 


0.0348(14) 


6 


0.7392(26) 


0.0316(19) 


7 


0.7707(27) 


0.0300(19) 


8 


0.8008(39) 


0.0304(24) 


9 


0.8311(37) 


0.0295(23) 


10 


0.8606(53) 


0.0212(29) 


11 


0.8819(60) 


0.0322(26) 


12 


0.9140(60) 





Table 1: The potential and force at (3 — 6.2. 



estimate of the asympotic value we used the 4 : 3 T-ratio to evaluate V(R). For R > 5, we 
find a small but statistically significant difference between the 3 : 2 and 4 : 3 T-ratio results. 
This we extrapolate to large T using the estimate of the first excited state energy referred to 
above. Statistical errors come from a bootstrap analysis of the variation over our 18 sample 
configurations, where the variational path combination is fixed by the 1 : T-ratio analysis. 
This avoids the possibility that the variational approach is influenced by the statistical error. 
Results are shown in Table [l]. 

We also fit the static potential for R>2 with a lattice Coulomb plus linear term: 

V(R) = C-^- + KR (8) 
Rl 

where 1 / Rl is the discrete Coulomb Green function. Using a fit to the previously determined 
values of V(R) and taking into account correlations in errors between potentials at different 
.R-values, we find the result of Table |2|. Here the error quoted is statistical; it includes the 
error from the data sample variation and the T-extrapolation described above, but not the 
systematic effect of using different fit functions. This value for the string tension enables us 
to set the scale by requiring a~ l \J~K = 0.44 GeV, so yielding a" 1 = 2.73(5) GeV. This implies 
that the potential has been measured to a physical distance, Ra, of 0.87 fermi. Using the 
conventional 2-loop perturbative relationship between A and a yields y/K/aA^ < 85.9(1.5), 
where the inequality arises because of the observed lack of asymptotic scaling at (3 = 6.2. 



These results are in agreement with a similar analysis on 20 lattices at the same (3- value [15 



16 . This confirms the conclusion of that work that finite-size effects are small for the 



E 


K 


0.274(6) 


0.0259(9) 



Table 2: Fit to the force for R > 1. 



potential on lattices with L > 20. A recent independent analysis on a 24 3 x 32 lattice flTl 



gives a result for the string tension which also agrees with ours. From a comparison of our 
result for the string tension with those at neighbouring (3 values, we find that asymptotic 
scaling is not obeyed. This can be understood because the lattice bare coupling is a poor 
expansion parameter Jl8| . The study of large lattices at larger (3 confirms this [|19] and shows 
that the use of a more physical coupling gives excellent agreement with perturbation theory 
to two loops, allowing A to be related to the string tension. 



2.2 Quark Propagators 

To calculate the quark propagator we need to solve equations of the form 

{A - kB)G{x, 0) = 



(9) 



where, for either action, A is purely local and B connects only nearest-neighbour sites. We 
have investigated the least residual (LR) and least norm (LN) variants |20| of the Conjugate 
Gradient (CG) algorithm, and an Over-Relaxed Minimal Residual (MR) algorithm jn|, [22| . 
We also implemented a simple preconditioning which involves left-multiplying Equation (|9]) 
by (1 + kBA~ 1 ). This decouples even and odd lattice sites so that Equation (j|) becomes 



(A - K 2 BA-'B)G cvcn 

Qodd 



^- 1 r/ odd + nA~ l BG even . 



(10) 
(11) 



We solve Equation ( |T0| ) for G even and substitute into Equation (|TT|) to obtain G odd . For all the 
algorithms and for both actions, we found that this red-black preconditioning gave a factor 
of between 2.5 and 3.0 decrease in the number of iterations required, without increasing 
significantly the number of floating-point operations per iteration. 



Writing 



M = A — k 2 BA~ 1 B 
( = r) eyen + kBA' 1 ^™ 



(12) 
(13) 



the MR algorithm is 



C-MG 



even 




repeat until convergence 

a 



(Mn,ri 



{Mr h Mn) 
a' = Qa (over-relaxation) 

^even _ ^-ycvcn , / 

n+i = n - a'Mri. 

We found that choosing the over-relaxation parameter, Q, to be 1.1 gave a gain from over- 
relaxation of a few per cent. This is in agreement with the data in which indicates 
that the gains from over-relaxation decrease as /3 increases. However, for both actions and 
throughout the physical regime investigated, MR typically converges to a solution in fewer 
iterations than CG, with each iteration taking about half the time. Therefore, we always 
use an over-relaxed MR algorithm with red-black preconditioning. 

In order to obtain an improved matrix element of an operator containing quark fields, in 
addition to using the SW action, the quark fields must be rotated according to @, [23| : 

q -► q' = (l - - 7 - 5) q (14) 

q -> q' = q (l + - T S) (15) 

where we have used the equation of motion for the quark fields to remove the bare mass 
from the rotation. The lattice covariant derivatives are defined as follows: 

5„ f(x) = \ (U,(x)f(x + A) - Ul(x - p)f(x - A)) (16) 



and 



/(*) D,= \ (f(x + jl)Ul{x) - f{x - jl)U„{x - fi)) . (17) 



2 

We can define the rotated SW propagator, G R (x,y), by p4] 



G R (x,y) = (1 - It D)G(x,y)(l + \ T D). 



This propagator includes the rotations (14) and (|15|) and so the evaluation of all correlation 
functions is performed in the same way as with the unimproved action; no further rotations of 
the operators are necessary. The evaluation of the rotated propagator, G R , is only marginally 
more complicated than that of G; solve the equation: 

(A - kB)G'(x, y) = VxiV + - Vx , y T D (19) 
then G R is obtained easily from G' by 



G R (x,y) = (l-^TD)G'(x,y). (20) 



Recently, there has been increasing evidence that lattice measurements can be enhanced 



through the use of non-local, or smeared operators [16|, |25|. The original attempts to 
use smeared operators for fermions are described in references [27 , |28| . For propagator 
calculations, we can use a source, rj, in Equation (||), which is extended over some non-zero 
spatial volume, called source smearing, or we can smear the solution, called sink smearing. 



Giisken et al. [25] have proposed a gauge- invariant smearing method based upon using the 



solution of the three-dimensional Klein-Gordon equation: 

(1 - K S D 2 )r] xfi = 6 xfi (21) 
for the source in Equation (^), where D 2 is the three-dimensional operator 

D 2 f(x) = £ (U 3 (x)f(x + j) + U}{x - j)f(x - ])) . (22) 

3=1 

Alternatively, sink smearing corresponds to solving 

(1 - k s D 2 )G sl = G R (23) 

where G R is the solution of Equation (^). The extent of the smearing is controlled by the 
single parameter, ks- 

We label propagators with a local source and sink, LL, those with a local source but smeared 
sink, SL, and those with a smeared source but local sink, LS. In practice, we smear at the 
sink, as this has the advantage that the smearing can be undone with a single multiplication 
by (1 — KsD 2 ). Since we can thus easily recover the LL propagator from the SL propagator, 
we only write out the SL propagator to disk, a very significant advantage for our particular 
machine architecture. Details of our investigations of the effects of smearing are reported 



elsewhere 29 



We have calculated LL propagators at k = 0.1510, 0.1520, 0.1523, 0.1526 and 0.1529 for 
the Wilson action, LL and SL propagators at k = 0.14144, 0.14226, 0.14244, 0.14262 and 
0.14280 for the SW action; the latter values were chosen to match roughly the pseudoscalar 
meson masses computed in the Wilson case. The boundary conditions used were periodic in 
space and antiperiodic in time. The scalar hopping parameter, ks, used in the computation 
of the SL SW propagators, was taken to be 0.180, corresponding to a smearing radius [f29| of 



approximately 2. The extra computational cost involved in using the SW action is difficult 
to quantify, being machine dependent. For the red-black preconditioned MR algorithm 
which we employed, each iteration took roughly 35% longer than for the Wilson action, but 
taking into account other fixed overheads, most notably I/O, the total elapsed time for a 
propagator calculation at a given pion mass was between 10% and 20% longer, whilst the 
memory requirement was some 10% more than in the Wilson case. 



2.3 Analysis of 2-Point Functions 



We obtain the amplitudes and masses of mesons (baryons) by correlated least-x 2 fits of a 
single cosh (exponential) function simultaneously to the appropriate forward and backward 
propagators. We require that the fitting region, [imin, t max ], satisfies the condition that, 
within the small range allowed by our statistics, changing t m ; n gives the same mass within 
errors, while t max is taken as large as the growing noise/signal ratio allows. We find that the 
intervals [12, 16] for local sinks and [9, 13] for smeared sinks are satisfactory in this regard. 
The correlated x 2 /dof varies between 0.3 and 4, which indicates both that we are taking 
correlations into account and that we are getting reasonable fits. 

The covariance matrix is estimated using a single-removal jackknife. For simple averages, 
such as the propagator, this is the same as the data covariance matrix. We calculate the 
jackknife covariance matrix for the timeslices used in the fit. We employ a singular value 
decomposition in the construction of our matrix inverses, to detect singular or nearly singular 
matrices. Because of our limited statistics, we use timeslices 12, 14, 16, 32, 34 and 36 for 
local sinks, and 9, 11, 13, 35, 37 and 39 for smeared sinks. We currently have too few 
configurations to deal with the k — k correlations directly in the fits. Therefore we do 
individual fits for each k value. We recover the effect of the k, correlations in our analysis of 
the errors. 

To estimate our errors, we use the bootstrap resampling method, a concise description of 
which has been given by Chu et al. [OT. From our set of 18 Monte Carlo samples of the 
probability distribution for the link variables, we construct a bootstrap sample by drawing 18 
configurations independently and with replacement. We perform exactly the same analysis 
on this sample of 18 configurations as on the original set, to obtain a bootstrap estimate of 
average quantities. We then build up the bootstrap distribution by drawing 1000 bootstrap 
samples with corresponding bootstrap estimates of averages. We reconstruct the covariance 
matrix for each bootstrap sample so as to take into account the uncertainty in the covariance 
matrix. We then bin the bootstrap determinations of the average to determine confidence 
limits. The quoted errors are obtained by requiring that the central 68% of the bootstrap 
values lie within the error bars. 

For every quantity fitted, we use the same sequence of bootstrap samples and store the best 
fit to the original set, together with the bootstrap determinations. Thus the i'th bootstrap 
estimate of any quantity, is the result of a fit to the i'th sample of 18 configurations. In this 
way we preserve information on correlations between physical quantities at, say, different k 
values. 



3 Hadron Spectrum 



3.1 Masses and Matrix Elements 

We present results for the pseudoscalar meson (P), vector meson (V), nucleon (N) and A 
using the following local interpolating fields: 

P = u^d (24) 

Vi = ujid (25) 

N = e abc (u a C l5 d b )u c (26) 

A„ = e abc {u a C lllU b )u c . (27) 

For the vector meson, we average our correlators over the three polarisation states, for the 
nucleon we average the 11 and 22 spinor indices of the correlator, and for the A we project 
out the spin-3/2 component and average over the four spin projections. 

The mass estimates in lattice units obtained for the Wilson action using local sources and 
sinks are given in Table || Also included are the amplitudes, A, and x 2 /dof obtained in the 
single cosh or exponential fits to the zero-momentum timeslice propagators: 

v-/w v,t, ^ i A h (e- mht + e- mh( - Lt -^), (mesons) , , 

5>0M* f ~ , 1 mht + T ,J , \ (28) 

x [ A h e mh \ t < Lt/2, (baryons) 

described above, where L t = 48 is the total time extent of our lattice. The corresponding 
results for the SW action using local sources with both local and smeared sinks are given in 
Table |. 

The Edinburgh plots for the LL data for both actions are given in ||. That for the SW SL 
results is plotted in Figure 0. The plots are broadly consistent, showing a trend towards 
the experimental value for mjqjm p with decreasing pseudoscalar meson mass. However, a 
comparison of the SW LL and SL plots reveals that the errors in the SL data are smaller 
than those of the LL data by at least 50%. This reduction is attributable to the better signal 
for the nucleon correlator using the SL data; comparing the results for the baryon masses 
at the lightest quark mass for LL and SL in Table ^, we note that the signal increases from 
approximately 2a to around 4cr from zero. For this reason, we employ SL data, wherever 
possible, in making our comparisons of the SW and Wilson hadron masses. 



3.2 Chiral Extrapolations 



We now consider chiral extrapolation of the hadron masses in order to estimate the lattice 
scale for the two actions. 
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Table 3: Hadron masses and amplitudes, in lattice units, for the Wilson action. The last 
column contains the values obtained by linear extrapolation to mp = 0. 
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X 2 /dof 
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1.4/4 


0.567 + 2 

0.078 + 2 
2.7/4 


0.549 + 2 

0.071 + 2 
6.3/4 


0.528 + 2 



Table 4: Hadron masses and amplitudes, in lattice units, for the SW action. The last column 
contains the values obtained by linear extrapolation to mp = 0. 
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Figure 2: Edinburgh plot for the SW action using sink smearing. 





a-^GeV) 


physical quantity 


\A/ 1 1 cnn 


O VV I ±J±J 1 


O VV IOJj 1 


m p 


2.77 + _£ 


2.63 + £ 


2.57 + 2 




2.39+;; 


2.50 + S 


2.45 + ^ 




2.48 + J 


2.40 + ;5 


2.33 +;j 


(m K * - m p )/(m 2 K - ml) 


2.0512 


2.06 + 2 


1.93 IS 




2.73 + I 

— 5 



Table 5: Scales determined from different physical quantities. 



From the bootstrap analysis, we find that there are strong correlations between the pseu- 
doscalar meson masses at different « values. Figure ^ shows correlated linear fits to the 
data for m 2 P versus 1/2k for both actions, at all five k values, confirming PCAC behaviour 
throughout the quark mass range used. The strong correlations presumably account for the 
remarkable linearity of our best estimates for m 2 P as functions of 1/2/c, given the size of the 
statistical errors. From the chiral extrapolation, we obtain 





= 0.15328 + ^ 


(Wilson), 




^crit 


= 0.14313 + 4 


(SW LL), 


(29) 


^crit 


= 0.14311+3 


(SW SL). 





The lattice scales obtained from correlated linear extrapolations of the other hadron masses 
to the chiral limit are presented in Table [|. For both actions, the scales derived from m p 
agree well with the scale from the string tension, whereas the scales from the baryon masses, 
whilst consistent with one another, are lower, by 2 — 3a for the SW data and more than 3a 
for the Wilson data. We have investigated quadratic chiral extrapolations, but find that we 
can attach no statistical significance to any difference between the extrapolations. 

The value of k corresponding to the strange quark may be estimated by assuming that the 
pseudoscalar meson mass obeys 

m 2 P (K 1 ,K 2 )=b P (^- + ^ —) (30) 

for valence quark masses corresponding to K\ and k 2 - We determine bp and K crit from the 
data for degenerate valence quarks, and match mp(K s , K cr it) to the physical kaon mass, taking 
the scale from m p . This gives 

k s = 0.1517 + 1 (Wilson), k s = 0.1418 + 2 (SW SL). (31) 
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Figure 3: m 2 P versus 1/2/c. 



3.3 Mass Splittings 



In an attempt to highlight any differences arising from use of the SW action, we present 
estimates for the vector-pseudoscalar meson and A-nucleon mass splittings, which in QCD- 
inspired quark models arise from spin interactions and therefore may be corrected at 0(a) by 
the spin term in the SW action. The suggestion that the vector-pseudoscalar mass splitting 
might be susceptible to lattice artefacts has been made previously by the APE collaboration 
in a comparative study of Wilson and staggered fermion actions, albeit at stronger 
coupling than the present work. 

Experimentally, for both light-light and heavy-light mesons, niy — nip is very nearly inde- 
pendent of quark mass: 0.57 GeV 2 (p - it), 0.55 GeV 2 (K* - K), 0.55 GeV 2 (D* - D), 
whereas for the charmonium system m 2 ji^ — m 2 c =0.72 GeV 2 , suggesting that for systems 
composed of equal-mass quarks, this quantity increases slowly with quark mass. Taking the 
scale from the string tension, the range of pseudoscalar meson masses for which we have 
data is 330 MeV to 800 MeV. Thus, it is of interest to study the quark-mass dependence of 
our data for m v — nip, with equal- mass quarks. In Figure f|, we plot the quantity m v — nip, 
calculated directly from the bootstrap masses, versus nip, for both actions. The Wilson data 
is consistent with previous work |32|] which indicated a negative slope, inconsistent with ex- 
periment at large quark mass. The larger statistical errors in the SW data leave open the 
possibility of a reduced dependence on quark mass, although it appears that there is still a 
tendency for the hyperfine splitting to decrease with increasing quark mass. We have shown 
in a related study that at heavier quark masses this is indeed the case; the SW estimate 
of the hyperfine splitting in charmonium is a factor of 1.8 larger than the Wilson estimate, 
although still smaller by a factor of roughly two than the experimental value. Using the 
string tension scale, the experimental value for light quarks corresponds to 0.075 in lattice 
units, consistent with the three lightest Wilson points and all the SW points. 

The A-nucleon mass splittings for the two actions, estimated from the bootstrap samples, 
are shown in Figure |5|. Using the scale set by the string tension, the experimental result 
of 300 MeV would translate to 0.11 in lattice units, a value broadly consistent with both 
the Wilson and the SW data. However, the errors in both data sets are too large for this 
quantity to discriminate between the two actions. 

It has been noted that the ratio (mx* — mp)/{ m2 K ~ m V) typically gives a value for the 



inverse lattice spacing which is lower than that obtained using other physical quantities [31 



Assuming that the vector meson mass is linear in the quark masses, with slope by and 
intercept ay: 

my(«i, « 2 ) = ay + — m P (Kx, « 2 ) (32) 
bp 
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Figure 4: my — m 2 P versus mf>, in lattice units. 
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we obtain the scales given in Table [| from 



m K , - m p _ m v (n s , K crit ) - m v (K crit , K crit ) _ by_ ^ 

TTlj^ Tfl^ mp^K s , Kcrit) m p{^criti Kcrit) bp 

We find that the scales for the two actions are in good agreement with each other, but our 
best determinations are more than one standard deviation below the baryon scales and more 
than two standard deviations below the string tension scale. 



4 Meson Decay Constants 
4.1 Pseudoscalar Decay Constants 

We determine the pseudoscalar decay constant, fp, through the matrix element of the fourth 
component of the lattice axial current: 

^(0|g(0) 7 475?(0)|P(p)) = fpmp (34) 

where our normalisation is such that the physical value of f n is 132 MeV. The factor Za is 
required to ensure that the lattice current obeys the correct current algebra in the continuum 
limit fl35fl. 



Because the signal for J2x(Ai(x, t)A\(0)) is unacceptably noisy, for both the Wilson and 
SW actions, we determine fp using the LL propagators by fitting the ratio 

Ex(A,(x,t)4(0)) gm| 

where mp and (0\P\P) are obtained from a separate least-x 2 fit to J2x(P{ x > t)P (fy)- The 
errors in fp are determined by a bootstrap analysis on the whole procedure. 

The measurements of fp reported in |J for the SW action displayed much larger errors 
than those for the standard Wilson action. For the SW action, we computed the correlator 
J2x(A4,(x, t)P^(O)), and thus we are also able to obtain f P by fitting the ratio 



5> 4 (*,t)pt (0 )> 

4=^ 1 ~ — — — — — — tanhmp(L t /2 — t), (36) 

52(P{x,t)P\0)) Z A (0\P\P) n 11 h K } 

x 

where mp and (0\P\P) are obtained from the fit to 52 x (P(x,t)Pi(Q)). 

We found that, for the light hadron sector using Wilson fermions, the use entirely of local 
operators provided a good determination of fp, and the introduction of SL and SS propa- 



gators offered no improvement p9 |. Therefore, we did not pursue the determination of f w 
using the SL propagators. 
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0.14 + 2 
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Table 6: Values of the pseudoscalar decay constant, in lattice units, and the ratio fp/m v . 
The last row in each table contains values obtained by a linear extrapolation to the chiral 
limit. 

In Table ^|, we present the values obtained for fp and for the dimensionless ratio fp/mv, 
using the Wilson action with the local axial current, and the SW action with the 'improved' 
axial current: 

q{x){l + - 7 - D)lM l ~ \l- D)q(x). (37) 

The final row contains results for each column after linear extrapolation in 1/k to the chiral 
limit. The measurement of fp for the SW action through the axial-pseudoscalar correlator 
is clearly much less noisy than that through the axial-axial correlator, and will be used in 
the following discussion. Our lattice results for fp/my vary only slowly with quark mass, in 
agreement with the experimental observation that f^/nip (0.17) is approximately the same 
as fK/ m K* (0.18). The chiral extrapolation of this quantity may therefore be more reliable 
than that of fp alone. 

In order to determine the physical values, the lattice results given in Table ^| need to be 
multiplied by the appropriate renormalisation constant, Z^ or Z A . 



Zj ~ 1 - 0.132/ 



(38) 



in perturbation theory. If we use the bare coupling constant as the expansion parameter, 
then ZY — 0.87 and we find, for the Wilson action, f w /m p = 0.18 ^ 3. However, reference [|Tj| 
proposes the use of an 'effective coupling', g eS , defined through 

£=^^(1 + + Odd,)). (39) 

where g is the bare coupling. At (3 = 6.2, g^ s ~ l-75#o yielding Zjf ~ 0.78, and we obtain 

The perturbative estimate |36| of the renormalisation constant Z^ is close to 1, 



Z C A ~ 1 - 0.0177# 2 ; (40) 

using the bare coupling leads to Z^ ~ 0.98, whereas the effective coupling gives Z^ ~ 0.97, 
yielding, for the SW action, f^/nip = 0.13 ^ \ in both cases. We note that the uncertainty in 
Z& due to the choice of the perturbative expansion parameter is about 10% with the Wilson 
action and only about 1% with the SW action. 

A recent non-perturbative estimate of Z^, based on the use of chiral Ward identities, gave 
the result Z% = 1.09(3), about 10% higher than the one-loop perturbative values quoted 
above |p7 |. This result was obtained from a simulation at (3 = 6.0 using a single value of the 



quark mass. Although it is expected that the dependence of Z^ on the lattice spacing and 
quark mass should be very mild, we feel that this expectation should be checked before the 
non-perturbative value is adopted in the present calculation. We note however that if the 
non-pertubative value of Z^ proves to be stable, the discrepancy between the result for f^/mp 
which we obtain using the SW action and the physical value is considerably reduced. This 
underlines the importance of reliable non-perturbative determinations of the renormalisation 
constants in order to get better estimates of the remaining lattice systematic errors, such as 
quenching. 

Figure ^| shows fp / my against mp in physical units for both actions, with the lattice spacing 
determined from m p , and Z^ computed using the effective coupling. Although the behaviour 
of this ratio with mp for the SW action is very encouraging, there is a clear discrepancy 
with the physical values in the overall normalisation. 

If we assume that the pseudoscalar decay constant obeys 

f p = a f + b f mp(K 1 ,K 2 ), (41) 
where m P (Ki, n 2 ) is given by Equation (p0|), then 

- 1 = —m P (K Sj K crit ). (42) 
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Figure 6: fpjmy against m 2 P , with lattice spacing determined from m p , and Za computed 
using the effective coupling. 
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Table 7: Values of the vector meson decay constant. The last row contains values obtained 
by a linear extrapolation to the chiral limit. 

We take a/ and bf from the fit to the data in Table and using k s values from Equation (|3T|) 
we obtain 

^-1 = 0.16 + 1 ° (Wilson), ^-1 = 0.25 + ^ (SW), (43) 

Jit Jtt 

compared to the experimental value of 0.22. Note that the ratio is less sensitive to uncer- 
tainties both in the renormalisation of the axial current and in the scale. 



4.2 Vector Meson Decay Constant 

The vector meson decay constant is defined by the relation 



m 



^(01^0)7^(0)11/) = -% (44) 

Jv 

where e M is the polarisation vector of the meson and Zy is the renormalisation constant for 
the lattice vector current [[3^]. We determine fy, using the LL propagators, by fitting to 

EE**)*) ~ —^e m ^ 2 coshm v (L t /2-t). (45) 

3=1 X ZZj vJv 

For the Wilson action we use the local current and for the SW action we use the improved 
local current. We obtain the values shown in Table 0. We extrapolate linearly in 1/k to 
obtain the values at the chiral limit. 

The renormalisation constants corresponding to our choice of the lattice vector currents are 
given in perturbation theory by 

Zy 7 ~ 1-0.17/ (46) 
Z$ ~ 1-0. 10g 2 . (47) 

The one-loop perturbative correction, although smaller for the SW action than for the Wil- 
son action, is still substantial and it introduces a significant uncertainty in f v due to the 
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Figure 7: l//y against m|>, with lattice spacing determined from m p , and computed 
using the effective coupling. 



uncertainty in the value of the expansion parameter. Recent non-perturbative estimates 



of Zy at f3 — 6.0 [37] give a value which agrees well with that obtained by using the ef- 



fective coupling defined in Equation (|39|). Since the non-perturbative values of Z v for the 
two actions are not yet known at /3 — 6.2, we use the effective coupling in the perturbative 
expressions for comparison with experiment, i.e., 



7 w 

J V 

7 c 

6y 



0.71 
0.83. 



(48) 
(49) 



We note that the perturbative uncertainty would be removed entirely by use of the conserved 
vector current. 

Our results for fy are compared with experimental values in Figure [7j. We do not know 
which pseudoscalar meson mass to associate with the experimental value of in this plot. 
For the purpose of illustration, we take the pseudoscalar meson mass to be the mass of the 
t]'. Despite this uncertainty and the presumably small residual uncertainty in the overall 
normalisation, we find the agreement with experiment encouraging. In the chiral limit we 



obtain 



l/f p = 0.33 1\ (Wilson) (50) 
l/f p = 0.31 (SW) (51) 

compared to the physical value of 1/ f p of 0.28(1). Again this discussion points to the need 
for a non-perturbative determination of the current renormalisation constants. 



5 Conclusions 



Significant differences have already been reported between the mass splittings obtained using 



the SW and Wilson fermion actions for systems involving heavy quarks [|33|, |3l|]. To obtain 
a complete picture of the effects of the 0(a)-improvement proposed by Sheikholeslami and 
Wohlert, it is necessary also to study the light quark sector. Our main conclusion is that at 
(3 = 6.2, for pseudoscalar meson masses in the range 330-800 MeV, there is no statistically 
significant difference between the results for the hadron spectrum or decay constants obtained 
with the two actions. This supports the validity of results obtained for light hadrons with 
the Wilson action over the last few years. 

We find that the scales obtained from the meson sector are consistent with that from the 
string tension, but there is evidence of inconsistent scales from the baryon sector and from 
the ratio (mx* — m p)/( m2 K ~~ m V)- 

Our calculation of fp/my with the SW action yields a result whose dependence on the quark 
mass is consistent with experiment, but whose magnitude is significantly smaller. There is 
a small residual uncertainty in the calculated values, which could be removed by a non- 
perturbative determination of the axial vector current renormalisation constant. It seems 
likely that this will not account entirely for the discrepancy and may signal the effect of 
quenching. Our results for f p and are broadly in agreement with experiment. 

We believe that this comparative study has established the viability of the SW formulation 
of quenched lattice QCD, with no significant disadvantages for light hadrons, laying the 
foundation for its use in the study of heavy-quark systems where discretisation errors are 
more serious. 
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